Quantum limited sensitivity of SET-based displacement detectors 
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We consider a model of a quantum-meclianical resonator capacitively coupled to a single electron 
transistor (SET) . The tunnel current in the SET is modulated by the vibrations of the resonator, and 
thus the system operates as a displacement detector. We analyze the effect of the back-action noise 
of charge fluctuations in the SET onto the dynamics of the resonator and evaluate the displacement 
sensitivity of the system. The relation between the "classical" and "quantum" parts of the SET 
charge noise and their effect on the measured system are also discussed. 



Micro-mechanical resonators have been used as ultra- 
sensitive force detectors in a number of experimental 
appUcations, ranging from Atomic Force Microscopy to 
Magnetic Resonance Force Microscopy [1] to experiments 
on Casimir force detection [2] . Recently, mechanical res- 
onators with vibrational eigenfrequencies (ly) of the order 
of 1 GHz have been fabricated [3]. At low temperatures 
{hv > ksT ~ 50 mK), these resonators provide an ex- 
ample of a man-made system that can be used to test 
the basic principles of quantum mechanics at the macro- 
scopic level. 

The standard cantilever displacement measurement 
schemes are based on laser interferometry, and can reach 
the levels of sensitivity of the order 10~'*A/a/Hz. This 
level of sensitivity requires, however, high laser power 
that may not be compatible with the ultra-low tempera- 
ture operation. This limitation provided the motivation 
to explore alternative electical measurement schemes 
[4-6]. In particular, Blencowe and Wyborne [4] have 
suggested based on a semiclassical analysis that by ca- 
pacitively coupling the cantilever to a Single Electron 
Transistor (SET), it is possible to achieve the sensitivity 
better than the zero-point-motion uncertainty. More re- 
cently, two of us [6] have found based on a fully quantum 
mechanical description of the quantum measurement of 
a cantilever using a quantum point contact (QPC), that 
the apparatus back action (current shot noise that in- 
duces force noise on the cantilever) fundamentally limits 
the displacement sensitivity and leads to a quantum-to- 
classical transition in the oscillator dynamics. Due to the 
resonant nature of transport through SET, it is expected 
to have significantly higher displacement sensitivity that 
QPC, and hence is more attractive from the experimental 
standpoint. Here, we analyze the fundamental sensitiv- 
ity limits of an SET-based detection scheme. We find 
that the higher "classical" sensitivity comes at the ex- 
pense of drastically increased back-action (island charge 
noise), which leads to large rms fluctuations of the can- 
tilever and to deterioration of the oscillator quality factor 
in the experimentally most attractive threshold regime. 
The optimal sensitivity is achieved in the co-tunneling 
regime in which the SET scheme becomes equivalent to 
the QPC-based detection. Similar conclusions have also 



been reached in the study of SET charge sensitivity [7] . 
We find that under no circumstances is it possible to ex- 
ceed the standard quantum limit using continuous SET- 
based detection. 




FIG. 1. Schematics for the model setup. 



The model is schematically presented in Fig. 1. The 
measuring apparatus is a single electron transistor: a 
quantum dot coupled to the leads via tunnel junctions. 
For simplicity we assume that the dot contains a single 
resonant level. The mechanical system, in the follow- 
ing referred to as an oscillator, which can be either a 
micro-mechanical resonator or a localized phonon mode, 
is capacitively coupled to the resonant level. Therefore, 
the displacement of the oscillator can alter the position 
of the resonant level with respect to the chemical poten- 
tials in the leads and thus can affect the tunnel current 
through the device. 

The Hamiltonian of the model can be written as 

H = Hloads + Hose + H' , (1) 

where the first two terms are the Hamiltonians of the elec- 
trons in the leads and the oscillator respectively, iJieads = 
'Eq.n=L,R.'^'lnc\r^Cqn and H^sc = im/2)dl + V (x) . Here 
Cq n=L(R) i'^q,n=L{R)) crcatcs (annihilates) an electron 
with a quantum number q in the left (right) lead, m and 
luq are the oscillator effective mass and frequency, and x 
is the oscillator coordinate. Here and in the following we 
set both Planck's constant h and the electron charge e 
equal to unity, unless stated otherwise. The single par- 
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tide states Eqn are filled up to chemical potentials in 
the leads, /xl and hr, which are biased by external volt- 
age, Hl — IJ'R = V, see Fig. 1. For simplicity we assume 
zero external temperature. The Hamiltonian H' includes 
both electron tunneling and modulation of the position 
of the resonant level by the oscillator: 



H' 



q,n=L,R 



(2) 



In the first, tunneling, term of the Hamiltonian (2), the 

operator d\ (d) creates (annihilates) an electron in the 
resonant level, ho = d^d, and the tunneling amplitude T„ 
is assumed to be independent of the single particle states 
in the leads. We assume that the energy of the resonant 
level Eo{x) depends linearly on the oscillator's coordi- 
nate X, i.e. Eo{x) = eo -|- Aa;. The unperturbed position 
of the resonant level eq can be varied by appropriately 
adjusting the gate voltage Vq, see Fig. 1. The parameter 
A physically represents an effective electric field in the 
capacitor formed by the oscillator and the quantum dot. 

We use the Keldysh-Feynman- Vernon formalism [8,9] 
to determine the evolution of the oscillator under the in- 
fluence of the tunneling electrons. We define a scattering 
operator for the oscillator alone, i.e. with electronic de- 
grees of freedom traced out: 

5osc = Troi [pei%S{- 00, oo)S{co,-oo)] /Tr [pci] ■ (3) 

In Eq. (3) S{oo, — oo) and »S(— oo, oo) are scattering oper- 
ators for the full system, »S(oo, — oo) = exp[— i J^^Hdt], 
where H is defined in Eqs. (1,2), and the operator % 
denotes time ordering along the Keldysh contour. The 
density matrix of the unperturbed electrons is the direct 
product of the uncoupled density matrices of electron 
reservoirs in the leads {pL and pu) with an empty electron 
state in the resonant level (po = dd^), pei = Pl^Pr^Pd- 
Eq. (3) implies that at t = — oo the leads, the resonant 
level, and the oscillator are uncoupled and that the in- 
teraction, H', is switched on adiabatically at t > — oo. 

In what follows we assume that the coupling constant 
A is small, while the tunnelling amplitudes Tl and Tr 
need not be small. Then, <Sosc can be written explicitly 
as a functional integral over the oscillator coordinate as 



«5osc = j 'C'Xc exp i dtC'^ 



X exp 
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dtidt2X{ti)x{t2)K{ti -t2) + 



■ (4) 



In this work we limit ourselves to O(A^) contribution to 
the effective action of the oscillator. The higher orders in 
the expansion (denoted by ... in Eq. (4)) are unimportant 
in the limit of strong tunneling as will be seen below. 

The first order contribution in A, i.e. interaction of 
the oscillator with average charge (no) in the dot, is 
included in the Lagrangian of the oscillator in Eq. (4), 



= £osc'^°'' — X{ho)cix. The charge (no) is related 
to the Fourier transform of the renormalized single par- 
ticle Green's function Gnit' - t) = -i{%d{t)d^t'))e\ 
as (no)ci = (l/27rz) J dujGj)^{u!). The averaging de- 
noted by ( )ei, is taken with respect to the exact sta- 
tionary state of the electronic subsystem alone, i.e. de- 
coupled from the oscillator. The rcnormalization of 
Goit) by the tunneling transitions can be obtained 
by a standard calculation [11]. Following the nota- 
tion of Ref. [8,10] we define a matrix Green's function 
G^(t2 — ti) = —i{d{t{)d^{t2))e\, where t\ and ^ can 
either be on the same or different Keldysh contours, 
i.e. i, j = ±. We also introduce the unperturbed 
Green's functions of the electrons in the left and the 
right leads G'^^^giti - t2) = -i{%c{i{)qnC^{t\)qn)e\,o, 
n = L,R, and the unperturbed Green's fimction of 
the dot electron Gg_o(^i " ^2) = -i{%d{t{)d'^ {t\))^ifi. 
The time ordered and anti-time ordered Green's func- 
tions, i.e. with time arguments on forward and re- 
turn branches respectively can be expressed in terms of 
the Green's functions with time arguments on different 
branches as G++{t) = e{t)G-+{t) + e{^t)G+-{t) and 
G—{t) = e{t)G-+{t) + e{-t)G+-{t), where e(i) is a 
unit step function [8,10]. Then, by solving the Dyson 
equation = q + G^qS'^'G^, where self-energy 

Y^ij _ J2q,n=L,R^n^qn,0' after straightforward calcula- 
tion we obtain: 

G^(a;)-2^ _ ^^^^ + (r^ + r^)^ ' ^^^^ 

= (^-eo)^ + (r, + r,)^ ' 

where G^(lj) = / G^^{t) exp {iu!t)dt. In Eq. (5) we intro- 
duced tunnelling rates ^l(r) = ''^T'^(j^)Pl{r)j where the 
densities of states in the leads Pl(r) are assumed constant 
for simplicity. 

The O(A^) contribution to the effective action in 
Eq. (4) is generated by the integral kernel K(t\ — ^2), 
which is just a two-point correlation function of charge 
fluctuations in the dot, and can be expressed as a product 
of two single particle Green's functions Go- The double 
integral in Eq. (4) can be rewritten as 

/ dtidt2[2ix''{ti)x'i{t2)Q{ti - t2)A{ti - t2) 
J —00 ^ 

+x<i{h)xi{t2)S{h-t2)}, (6) 



where we have introduced the "rotated" Keldysh vari- 
ables x'^it) = x{t+) - x{t-), x''{t) = xif+) + x{t-). The 
kernels A(ti - ^2) = Im[G^"(t2 - ti)G^+(ti - ^2)] and 
^(ti-ia) =Re[G+-(<2-<i)G^+(ti-t2)], with G+- and 
G^)^ given by Eqs. (4), are related to anti-symmetric 
(quantum) and symmetric (classical) parts the charge 
correlation function K{t2 — ti). 
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If we assume that the tunneling through the resonant 
level is fast compared to the motion of the oscillator, i.e., 
the kernel K{t2 —h) is nonzero on a time scale which is 
much smaller than the time scale of the oscillator, then 
the oscillator trajectory between ti and t2 is close to a 
straight line, i. e., x''^i\t2) =i x''^i\ti)+x<i\ti){ti-t2). 
Using this approximation in Eq. (6) gives an action which 
is local in time: 



\2 n 

T=-— dt {iRx'^xi + iAx'^xi + Sx'ix") + . 



(7) 



where the the coefficients R, A and S can be expressed 
in terms of single particle Green's functions as follows: 



S=^J duG-^+{co)G+-{u), 

G^+(a;i)G+-(a;2) 



R= —V I dioidLo2 

TT 



UJl — (jJ2 



(8a) 
(8b) 
(8c) 



The effective action in Eq. (7) is exactly of the form 
of the Caldeira and Leggett action [12], derived for a 
bosonic heat bath at high temperature. The first term 
in the effective action is a rcnormalization of the 
oscillator potential. In contrast to the infinite rcnor- 
malization in the Caldeira-Leggett model, here R is fi- 
nite. Evaluating the integral in Eq. (8c) using Eqs. (5) 
yields R — 2d{no)ei/deo, where (no)oi is the average oc- 
cupation of the electron in the dot given by Eq. (6). 
Physically this corresponds to the back action of tun- 
nel current, which is perturbed by the displacement of 
the oscillator. Combining this second order rcnormal- 
ization with the first order rcnormalization in Eq. (4) 
the effective potential of the oscillator can be written as 
V{x) ~ V^^^'^ix) + Xx{na{x))ei, where (no(a;))oi is the 
occiipation number of the resonant level for a fixed po- 
sition of the oscillator. The last term in the exponent 
of Eq. (7) provides dephasing of the system. Its effect 
at the classical level corresponds to a white noise force 
f{t) exerted by the tunnel current on the oscillator. The 
second term causes energy damping. The classical equa- 
tion of motion for the oscillator [12] can be written as 
mx + m'yx + dxV — f{t), where, in our case, m7 = X^A 
and if (t) fit')) = X'^SSit - t'). Thus the classical and 
quantum parts of the resonant level charge correlation 
function determine fluctuations and dissipation for the 
oscillator respectively. One can therefore define an ef- 
fective temperature, Tee, using a fluctuation-dissipation 
relation, giving S/A = 2Tes- The effective temperature, 
Tefi, is not determined by the reservoir's actual temper- 
ature, as in the Caldeira-Leggett model [12], but rather 
by the coupling to the tunnel current. 

Tefi' determines the fluctuations of the oscillator coordi- 
nate due to the tunnel current induced noise. In the case 
of a linear oscillator, V{x) = muil^^x^ /2, the dispersion 



of the oscillator coordinate is {x"^) = T^a / {itiluI^^) . We 
can now check the validity of our expansion in A. From 
the structure of Eqs. (7,8), we see that higher order terms 
in Eq. (4) will be smaller by powers of the dimensionless 
parameter X^Jix^ /{Tl + ^r)- Physically, if the oscilla- 
tor induced shift of the resonant level is small compared 
to the width of the level, the back action on the oscillator 
is weakly dependent on the position of the oscillator and 
the higher order nonlinearities are unimportant. Thus, 
for sufficiently large F, we only need to consider the lead- 
ing, quadratic, terms in the effective action (7). 

It is instructive to evaluate Teg explicitly, using 
Eqs. (5,8) in two limiting cases-the threshold and co- 
tunneling regimes. Suppose first that the resonant level 
is in the vicinity of one of the chemical potentials in the 
leads, say eo = = Q and the bias between the chemical 
potentials is large, /il = oo. In this regime, the current 
through the device is very sensitive to the energy of the 
resonant level, and hence to the displacement of the can- 
tilever. In this case we obtain: 7'^'' = hX^T J {■KmV^) ^ 
Tlfi = 7rFi/4. Therefore, the effective temperature of 
the oscillator in the threshold regime is essentially de- 
fined by the tunneling-induced width F of the resonant 
level. For a practical Si nano-mechanical resonator with 
dimensions 3/im x O.l/xm x O.l/xm, coupled to SET with 
F ^ 10^ s~^ by a coupling strength A ^ 10~^V/A, this 
corresponds to an effective temperature T^s ~ 0.1 K, and 
the damping coefficient 7 ~ 10^ s~^. The effect of back 
action in this case will limit the lowest achievable oscil- 
lator temperature to 0.1 K, and the maximum quality 
factor to about 100. 

Another important limiting regime corresponds to the 
situation when the resonant level is far above or below 
the chemical potentials in the leads. This is the so-called 
co-tunneling regime. The tunneling electrons can now 
occupy the level only virtually and the effective coupling 
constant between the leads is small as it is suppressed 
by the large energy separation between the chemical po- 
tentials in the leads and the resonant level. Assuming 
Ml - Mfl = ^ < (ml + I^b)/2 = fiF, and i^l{R) » P, eo, 
we obtain 7™* = hX'^T^ / {^Tm^4) , r<f|* = TlTrV/T'^. 
Compared to the threshold regime, in the co-tunneling 
both the oscillator damping and the back-action noise 
are significantly reduced, with the bias voltage across the 
leads determining the effective temperature of the oscil- 
lator. This result is consistent with Ref. [6,13]. 

We are now in a position to analyze the sensitivity of 
the system. Suppose that the oscillator (which will be 
assumed linear from now on) is perturbed by an external 
force F{t), say, a short kick of duration tf ^ uJq^ , so that 
F{t) ~ FrpSit). This kick results in the variation of the 
oscillator's amplitude by the amount Sx = Ftf / (jnuio). 
What minimum 5x can be detected by observing the tun- 
nel current, given the noise ([/u]^) in the current? 
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FIG. 2. Sensitivity in dimensionless units {xq/ 5x^^^) as 
a function of chemical potentials for a symmetric structure 
{Tl = = r/2). The axes are ^l/V and ^r/V. 

The ability to measure a signal can be represented by 
the integrated signal-to-noise ratio [14] 

5/n=(l/27r) j M{SF{i^))?im'). (9) 

where S'f(w) is the Fourier transform of the detector's re- 
sponse to an external perturbation, i.e., the force F. The 
variation of the current through the structure due to the 
variation of the oscillator coordinate is 51 — {dl /dx)Sx, 
and therefore the response can be written as Sf{uj) — 
{dI/dx)F{uj)/[m{Lu^ + i-foj - oj^)]. 

The current and the noise can be easily evaluated 
if we recall that the dynamics of the oscillator is 
slow compared to the dynamics of the tunneling elec- 
trons. The current in the adiabatic approximation, 
i.e., for a fixed position of the oscillator is given by 
I{x) = (l/27r) J dtoT {lo , x) ^ where the transmission 
coefficient T{uj,x) = ATlTr/[{uj - Eoix))"^ + F^] [11]. 
The noise at low frequencies (of order ujq) is given by 
~ (|/o|2)-^(a//c'a;)2(|Aa;^|2), where the shot noise 
is related to the transmission coefficient as (|/o|^) = 
(l/2^)/;;da;T(u;,0)[l-r(u;,0)] [15], and {\AxJ') = 

2'yTcff / m[{uj'^ — ^o)^ + 7^"^^] is the fluctuation spectrum 
for the oscillator. Substituting these expressions into 
Eq. (9) and setting the s/n = 1 as a criterion for a suc- 
cessful measurement, one obtains a criterion for a detec- 
tion of a minimum force by our apparatus. By expressing 
this minimum force in terms of the minimum displace- 
ment 5xmin that it causes, this criterion reduces to 

xl/6xl,, ^ \dI/deo\/ {^^/S(M)) , (10) 

where Xq — h/{2mujQ) is zero point displacement for the 
oscillator, and S is given by Eq. (8b). 

The sensitivity defined by Eq. (10) can be evaluated 
and is presented in Fig. 2 as a function of chemical po- 
tentials in the leads relative to the position of the reso- 
nant level. The sensitivity is maximal in the co-tunneling 
regime, where it reaches 1/2. In the threshold regime 



the sensitivity is somewhat smaller (~ by a factor of 
2). The sensitivity is worst when the resonant level is 
positioned symmetrically with respect to the chemical 
potentials (the blue regions in Fig. 2), as the current 
sensitivity dl/dx vanishes in this regime. These results 
are similar to the conclusions reached in the studies of 
the "non- ideality" of SET detectors [7]. Qualitatively, 
the reduction of sensitivity in the sequential tunneling 
regime can be attributed to the detector latency during 
the electron dwell time on the island, which contributes 
to the back-action noise, but not to the measurement. 

In summary, we have analyzed the quantum measure- 
ment of a mechanical oscillator coupled to an electronic 
resonant level that models a single electron transistor. 
We determined the back action effects of the detector 
on the quantum system, which lead to a measurement- 
induced effective temperature and damping coefficient. 
We also determined the fundamental sensitivity limits of 
the scheme in all operation regimes. 
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